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Abstract 

The Reynolds number dependence of the dimensionless dissipation 
rate C £ = eL/U 3 is derived directly from the Karman-Howarth equa- 
tion as C e = C e>0 o[l + A/Rl], where the coefficient A depends on 
the second- and third-order structure functions, Rl = UL/u, L is the 
integral length scale, and U is the rms velocity. Fitting this form to 
results from DNS indicates that A is effectively constant and hence 
the predicted dependence of C e on Reynolds number is as R^ 1 . 

PACS 47.27.Ak, 47.27.E- 

In recent years, there has been great interest in the Reynolds number 
dependence of the dissipation rate in homogeneous, isotropic fluid turbulence: 
see [TJ- [TTJ . Apart from its intrinsic fundamental interest, this is a key factor 
in the free decay of turbulence, which is one of the most studied aspects of 
the turbulence problem: see, for example, [12]- [15], and the many references 
therein. 

We consider the dimensionless dissipation rate 

e = C £ U 3 /L, (1) 

which was put forward in 1935 by Taylor [TB] on the basis of dimensional 
arguments. Here U is the rms velocity of the fluid and L is the integral length 
scale. As early as 1953, Batchelor [17] (in the first edition of this book) 
presented evidence to suggest that C e tends to a constant with increasing 
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Reynolds number. However, the present interest in the subject stems from 
the seminal paper by Sreenivasan [I], who established that in grid turbulence 
C e became constant for Taylor- Reynolds numbers greater than about 50 but 
observed that the actual value could depend on the flow configuration and 
initial conditions. 

Attempts have been made to establish a theoretical relationship between 
the dissipation rate and the Reynolds number. Lohse [32] used a mean-field 
closure of the Karman-Howarth equation to obtain an approximate form, 
whereas Doering and Foias [19] have established upper and lower bounds 
to be satisfied by such relationships. In this Letter, we derive a new and 
exact relationship between the dissipation and the Reynolds number. We 
begin by stating the Karman-Howarth equation and reviewing the relevant 
phenomenology. 

As is well known, the Karman-Howarth equation [20] is derived directly 
from the Navier- Stokes equation and is an exact relationship expressing con- 
servation of energy. It may be written in terms of structure functions as 
[21]: 

J_d_ ( v_b_( .dSA 

3 £ 2 dt 6r 4 9r lr 3j r 4 dr \ dr J ' U 

where the structure function of order n is given by 

S n = <(«(r) - «(0)) B > . (3) 

We begin by re-writing it as an expression for the dissipation rate, thus: 

= _3 0S2 _ J_ d_, 4 , to_d_ ( .dSA 

4 dt 4r 4 dr l? ~ 3> 2r 4 dr \ dr J ' U 

Now we re- visit the Richardson- Kolmogorov phenomenology of this equa- 
tion [T7], as this will be of specific use to us later on. We first consider the 
stationary case and, as usual, this requires the injection of energy at scales 
r > rj (say). As the Karman-Howarth equation is local in r, we restrict our 
attention to scales less than the injection scale rv and set the time- derivative 
equal to zero. With these steps, equation fll]) becomes: 



J__9 

4r 4 dr 




(5) 



We note that each term on the right hand side is separately a function of r 
but that jointly they are constant, as required to match the dissipation rate 
on the left hand side. If we go further, and restrict the values of the scale 
to r > ra, where is the dissipation lengthscale, then for sufficiently large 
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values of the Reynolds number the well known arguments of Kolmogorov 
(K41) [221 123] apply- In order to apply this rigorously, we take the infinite 
Reynolds number limit. This means taking the limit v — » such that e 
remains constant [17] . The viscous term may then be set equal to zero, 
and the remaining term on the right hand side represents the inertial flux of 
energy, which is now constant, and equal to the dissipation rate. Hence: 

From this equation it is a simple matter to recover Kolmogorov's famous 
'4/5' law for the third-order structure function, thus: 

S 3 ( r ) = —sr. (7) 
5 

Next we make a change of variables based on constant velocity- and 
length-scales, V and /, respectively. As an example, we will put the second- 
order structure function into a dimensionless form. We do this by expanding 
5*2 in a Taylor-Maclaurin series. Note that this is a general operation and is 
not restricted to any particular values of r. Then, introducing a dimensionless 
variable x, through r = Ix, we have 

S 2 ( r ) = S 2 (0) + S' 2 (0)r + ^ 2 '(0)r 2 + . . . , 

= 5 2 (0) + S' 2 (0)lx + ^S' 2 \0)l 2 x 2 + . . . , 



V 2 

-2 



/ 2 (0) + ^(0)x + ^'(0)x 2 + 



= V'Mx), (8) 

where, in the third line, we have renamed coefficients in order to extract 
a factor V 2 , and also to absorb powers of I. Note that the primes denote 
differentiation with respect to r in the first two lines but with respect to x in 
the third. Lastly, we have defined the dimensionless second-order structure 
function f 2 (x) as the sum of the Maclaurin series in the new coefficients. 

We can repeat this process for structure functions of any order n, and we 
summarise the general result as: 

S n = V n f n (x) } with x= r -. (9) 

We should note that this procedure involves no approximations or non-trivial 
assumptions. It merely introduces the f n as dimensionless forms of the struc- 
ture functions. As they are dimensionless, their dependence on r must be 
scaled by some length, here denoted by /. 
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Then we make the specific choices V = U, the root-mean-square velocity, 
and I = L, the integral length scale. With these choices, and substituting 
from 09]), we change equation (EJ) to the form: 



A 3 U 3 A 2 uU 2 



+ 



L 2 



where the coefficients A 3 and A2 are given by 

1 d 



A, 



4x 4 dx 



x f 3 (x) 



and 



Ao = 



_3_d_ 

2x A dx 



x 



dx 



Then, if we divide both sides of ( TTUj) by U 3 /L we may write this as 



C £ = A 3 + 



Rj 



A, 



A, 



A 3 R L 



(10) 



(12) 



(13) 



where the dimensionless dissipation C £ is as defined in ([I]) and the Reynolds 
number is given by Rl = UL/u. Note that the are determined by 
along with this choice of scaling variables, and hence so also are the coeffi- 
cients A 3 and A2. We also note that this equation is still just the Karman- 
Howarth equation: no approximation has been made. 

Let us consider the asymptotic behaviour of this expression for large 
Reynolds numbers. In general the coefficients A 2 and A 3 may depend on 
the Reynolds number (although, as we shall see later, comparison with the 
results from DNS suggests that they are essentially constant, albeit with 
possibly some dependence on initial conditions). However, the asymptotic 
properties of ([TBI must be those of the Karman-Howarth equation. So taking 
the limit of infinite Reynolds numbers, and comparing fflBl with ([5]), we note 
that the coefficient A 3 becomes constant, while the term A 2 /Rl vanishes. 
Then we obtain, by analogy with ([H]), the asymptotic form of ffT3"j) as 



lim C £ = C £ 
Substituting back into (TIB"]) then yields 



A, 



where 



C £ = C £:OQ [1 + A/ Rl] 



A = A 2 /A 3 



(14) 



(15) 



(16) 
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We now discuss the extension of our results to freely decaying isotropic 
turbulence. Of course, the neglect of the time-derivative term in ([2]) is quite 
usual, even for decaying turbulence, provided that the Reynolds number 
is large enough and that one restricts attention to the inertial range. For 
instance, this step is required in order to derive the '4/5' law for decaying 
turbulence and is known as local stationarity. However, we wish to consider 
a more general approach in which we introduce the time dependence of all 
statistical quantities, so that we now have e(t) and S n (r,t), at any time t. 

In order to make comparisons with the stationary case, we take some 
fiducial time t = t e , when the turbulence is assumed to have evolved from 
arbitrary initial conditions to a state determined solely by the Navier-Stokes 
equations. Then, as before, we introduce a change of variables, such that 

S n (r,t) = U n (t e )g n (x,r), (17) 

where now 

x = r/L(t e ); r — t/T; T — L(t e )/U(t e ). (18) 

Just as before, in the stationary case, the introduction of the dimensionless 
(but now time-dependent) structure functions g n (x,r) may be accomplished 
by the use of Taylor series; although, in this case, it is for a function of two 
variables. 

It should be emphasised that the form ffTTl) is not a similarity solution of 
any kind. If we were to drop the dependence on the variable r then this would 
amount to an assumption of self-preservation, as introduced by von Karman 
in 1938 [20J. However, this is one of the vexed issues of turbulence theory. 
It runs into difficulties, not least because during free decay the characteristic 
length changes from the integral length scale to the viscous scale in the final 
period of the decay: see, for instance, [13]. As it is, we reiterate that we do 
not make any such assumption and that we retain the full time dependence 
of the problem. 

With this in mind, it is easily shown that substituting (j!7p into (jl]) leads 
to a generalization of equation ([TBI to the form: 



c, = (A, - b 2 ) 



i + 1 A ' 



Rt — Z?9 



(19) 



where the coefficients A 2 and A 3 are still defined by equations ( TT2|) and (TTT|) . 
but now with / replaced by g, and the new coefficient B 2 is given by 
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Let us now consider the effect of including this time dependence. We 
introduce the inertial flux which we denote by H max - (Formally, this quantity 
is defined in wavenumber space and for our present purposes we are interested 
in its maximum value with respect to wavenumber: for a discussion see |24j.) 
Then, as is well known, for increasing Reynolds number, the maximum flux 
approaches the dissipation rate from below; or: 

lim U max ->• e, (21) 

Rl— >od 

for stationary turbulence. Our present analysis indicates that this cannot, 
in principle, be the case for freely decaying turbulence. 

We may see this as follows. Rewriting the coefficient B 2 in terms of the 
structure function, we have 



n 3 002 

B2= TWF 



3 L(t e ) dS 2 



(22) 



AU 3 (t e ) dt 

For free decay, this coefficient must be negative. Taking its modulus we have 

lim e^^-(A 3 + \B 2 \). (23) 
Thus, in freely decaying turbulence, we have 

lim > r—r < 1, (24) 

Rl->oo e A 3 + \B 2 \ 

and so the inertial flux of energy never quite reaches the same value as the 
dissipation. This result has implications for the interpretation of the Kol- 
mogorov (K41) picture, when a comparison is made of forced and decaying 
homogeneous, isotropic turbulence, and we shall investigate this further in 
future work. 

For completeness, we now give the extension of equation (TT3|) to the case 
of free decay, as: 

C £ = C £i0O [1 + A/R L ] , (25) 

where the dimensionless dissipation at infinite Reynolds number and the 
coefficient A are now given by 

C £>00 = A 3 - B 2 and A = A 2 /(A 3 - B 2 ). (26) 

We may compare our results to those of other theories. Lohse [18] used 
a mean-field approximation to the Karman-Howarth equation and obtained 
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the asymptotic result C £jOQ = (6/6) 3//2 , where b is the prefactor in the Kol- 
mogorov inertial-range form of the second-order structure function. This 
may be compared with our equation (JTlj) . which in contrast gives C £t00 = A3. 
This difference from our result is only to be expected because a closure in- 
variably involves expressing S3 in terms of S2. (It should be emphasised that 
our remarks here have no implications for the validity or accuracy of Lohse's 
approximation.) 

It is also of interest to compare equation (jlOp . which is an intermediate 
stage in our calculation, to the result for an upper bound on the dissipation 
as given by Doering and Foias [19]. This is featured in their abstract as 

U 2 U 3 
e <c 1 u— + c 2 —, 

and corresponds to their equation (40). Here the coefficients c\ and C2 de- 
pend on the shape of the forcing function, while / is its longest length-scale. 
Obviously this is quite different from our own result, where the correspond- 
ing parameters depend on the fluid turbulence and not on the forcing. Also, 
we have an equality, rather than an inequality. Nevertheless, it can easily be 
shown that the two results are equivalent and this will be presented later, in 
a fuller account of this work. 
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Figure 1: Values of C £ from DNS for two different initial spectra plotted 
against the integral length scale Reynolds number, TZl- These points were 
fitted to the curve y(x) — a + (b/x), as given by Eq. (l25]h 

Lastly, we illustrate the dependence of the dissipation on Reynolds num- 
ber by fitting the general form of equation fl25|) to the results of direct numer- 
ical simulations (DNS) of the Navier-Stokes equation. This is shown for two 



7 



Initial spectra 


Ci 


c 2 


c 3 


c 4 


Spectrum A 


0.0017 


4 


0.08 


2 


Spectrum B 


0.08 


2 


0.0824 


2 



Table 1: Initial energy spectra parameter values for use in the numerical 
computations. These parameters are substituted into Eq. f[2"T|) to generate 
the required initial spectrum. 

different initial spectra, which are generated using the following equation: 

E(k, 0) = C x k C2 exp { - C 3 k° 4 }, (27) 

where the values taken for the parameters C1-C4 are as given in Table [TJ Note 
that Spectrum A and Spectrum B have low-wavenumber regions varying 
as fc 4 and k 2 , respectively. The fitting process gives asymptotic values of 
C £j00 = 0.43 ± 0.01 and C £i00 = 0.43 ± 0.02, respectively. 

Figure [1] illustrates the way in which our new exact relationship between 
the dimensionless dissipation rate and the Reynolds number can be used 
to fit the data generated by numerical or other experiments. In particular 
it supports the idea that the coefficient A is constant with respect to the 
Reynolds number, even at low Reynolds numbers. 

Summing up, we have taken the Karman-Howarth equation as an ex- 
pression for the dissipation rate and transformed it into a more useful form, 
by making a change of variables. Then we have used its known properties 
to identifiy the asymptotic dimensionless dissipation coefficient and put it 
into the yet more useful form of fTTS]) . At no time have we made any ap- 
proximation or similarity assumption. Our only mathematical assumption is 
the underlying one of theoretical physics: that all variables corresponding to 
physical observables are mathematically well-behaved. From our comparison 
with DNS, we conclude that the low- Reynolds number dependence of the 
dissipation coefficient is of the form: constant x -Re- 
currently we are working to extend the simulations to forced turbulence, 
in order to study the stationary case. In particular, we intend to study the 
detailed behaviour of the coefficients at small Reynolds numbers, including 
the extent of their dependence on the initial conditions. This work will be 
the subject of a full report in due course. 
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